Constraints on Dark Energy Models Including Gamma Ray Bursts 
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In this paper we analyze the constraints on the property of dark energy from cosmological ob- 
servations. Together with SNe la Gold sample, WMAP, SDSS and 2dFGRS data, we include 69 
long Gamma- Ray Bursts (GRBs) data in our study and perform global fitting using Markov Chain 
Monte Carlo (MCMC) technique. Dark energy perturbations are explicitly considered. We pay 
particular attention to the time evolution of the equation of state of dark energy parameterized as 
wde — wo + w a (l — a) with a the scale factor of the universe, emphasizing the complementarity of 
high redshift GRBs to other cosmological probes. It is found that the constraints on dark energy 
become stringent by taking into account high redshift GRBs, especially for w a , which delineates the 
' evolution of dark energy. 
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Cosmological observations, including Type la supernovae (SNe Ia)[l], 0, d, Cosmic Microwave Background 
£SJ ■ radiation (CMB)|| @, Large-Scale Structures (LSS)0j and so on, have provided strong evidence for a spatially 
t-H | flat universe being in a stage of accelerating expansion. In the context of Friedmann-Robertson- Walker cosmology, 
this acceleration is attributed to a new form of energy with negative pressure, dubbed dark energy (DE), whose 
nature remains a big puzzle. The simplest candidate for DE is the cosmological constant. However, it suffers from 
the notorious fine tuning and coincidence problems [1, [§]. Many dynamical models on DE, such as quintessence 
[IS EH, E3, El , phantom [3 , k-essence fTsl Il6l] and quintom, have been proposed to avoid above difficulties. Among 
them, the quintom model allows the equation of state (EOS) crossing -1 durin g ev olution 11711 . Extensive studies on 
(N) ■ this model has been carried on both theoretically and phenomenologically [13, 111 GJ, H (H, SI H, HE HE S3, 

; [11, ii, [13, [M 111 [M HHi, m, H3, HI. 

Different DE models predict different global evolutions as well as different structure formations of the universe. 
Therefore cosmological observations can provide important constraints on the nature of DE. However, degeneracies of 
cosmological parameters exist in almost all cosmological observations, i.e., they are sensitive not to single parameters 
but to some specific combinations of them. It is therefore highly necessary to combine different probes to break 
parameter degeneracies so as to achieve tight cosmological constraints. Furthermore, different observations are affected 
' by different systematic errors, and it is thus helpful to reduce potential biases by combining different probes. To 
perform joint analyzes, global fitting is the most secure study, because it avoids using strong, and sometimes even 
d • inappropriate priors obtained from other observations. 

Observations on SNe la have played important roles in DE studies]!, [H, However, due to the difficulties of 
detecting high-redshift SNe la, such a geometrical measurement is limited. The future SNAP/JDEM is planed to 
probe up to redshift around z ~ 2. On the other hand, the 3-year WMAP data (WMAP3), which is so far the 
most precise measurement of the CMB, reveals the information obout our universe at redshift around z ~ 1100. In 
between, there have not been many probes accessible to us. In this regard, GRBs are the most promising tracers of 
the evolution of the universe at redshift around a few to even a few tens because they are the most powerful events 
in the universe. The currently operating Swift satellite is able to detect about 100 long duration GRB events within 
one year, and we expect that it will open up a potentially new window for GRB cosmographic studies. At present, 
the measured redshifts of GRBs have extended to z — 6.29 39]. Even though the physical origin of the long GRB is 
not very clear, many extensive discussions about the relations between the spectral and temporal parameters show 
the potential for using long GRBs as cosmic candles for cosmography. 

Recently, headway has been made in considering how to make GRBs standard candles. In the literature there are 
many studies about the intrinsic correlations between temporal or spectral properties of GRBs and their isotropic 
energies and luminositiy, for example, the spectral lag — the luminosity correlation (r — £i SO )[13], the luminosity 

— variablity correlation (V-Li SO ) 41], the spectral peak energy — isotropic energy correlation {E pea k-Ei SO ) 42], the 
peak spectral energy— isotropic luminosity correlation (Eneak-Li SO )^^ and the isotropic luminosity — peak energy 

— high signal time scale correlation {Li SO -E pea k-TQAb)uM an d so on - These correlations help GRBs nearly to be 
"known candles" [45[. Shaefer obtained the first GRB Hubble diagram of 9 GRBs with known redshift by using 
the spectral lag and the variability indicators, and from the GRB hubble diagram he constrained Q m < 0.35 within 
la confidence leve l 1461 . After that, many investigations have been triggered on using GRBs as cosmological probes 
[H [H [H, H3, HH [Ha [H, [H, HI . Very recently, Shaefer[13 has constructed a GRB Hubble Diagram with 69 GRBs 
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over a redshift range of 0.17 to > 6, with 39 GRBs having redshifts z > 1.5. It is the largest GRB sample so far. 
The aim of this paper is to present a more general analysis on the EOS of DE by including GRB data in addition to 
CMB, LSS and SNe la in global fitting. We mainly consider the dynamical dark energy parameterization for there 
is no compelling reason to neglect the evolution of dark energy factitiously. We employ MCMC techniques for the 
analysis. The MCMC method is the global fitting on the cosmological parameters and we use the original CMB 
and LSS data rather than use the CMB shift parameter, the linear growth factor or the Baryon Acoustic Oscillation 
(BAO) measured from the Large scale structure survey. The global fitting is the most secure way for processing data 
because it is a joint analysis and it can avoid using some strong or even inappropriate priors from others. The effects 
of dark energy perturbations are carefully taken into account with great attention paid to the perturbations when 
the equation of state gets across -1. Our paper is structured as follows: in Section II we describe the method and the 
data; in Section III we present our results on the determination of cosmological parameters from global fitting; finally 
we present our conclusions in Section V. 



II. METHOD AND DATA 



In this section we describe the method used in the fitting process. For the DE parametrization, we adopt KCDM 
model, constant w and pay particular attention to the commonly used EOS of the form[57l]: 

WDE(a) = w + w a (l - a) (1) 

where a = 1/(1 + z) is the scale factor and w a characterizes the "running" of the EOS. The evolution of DE density 
is then 

px(a)/px(a a = 1) = a -^ 1+wo+w ^ exp [3w a (a - 1)] . (2) 

For DE models whose EOS is not equal to —1, the perturbation inevitably exists. The perturbation of DE has no 
effect on the geometric constraints of SN la, however, for the CMB and LSS data, the perturbation of DE should 
be considered, because the late time ISW effects differ significantly when DE perturbation are considered, and the 
ISW effects take an important part on large angular scales of CMB and the matter power spectrum of LSS 22]. For 
quintessence-like or phantom- like models where w does not cross —1, there are not fundamental difficulties in dealing 
with perturbations. In parameter fitting, however, biases may be introduced if we limit our considerations only to 
quintessence or only to phantom models. Thus it is more natural and consistent to allow w crossing — 1 in the fitting 
analysis. Furthermore, there are observational indications that w might evolve from w > — 1 in the past to w < —1 at 
present, which have stimulated many theoretical studies. For w crossing — 1, one is encountered with the divergence 
problem for perturbations at w = — 1. For handling the parametrization of the EOS getting across -1, firstly we 
introduce a small positive constant e to divide the full range of the allowed value of the EOS w into three parts: 1) 
w > — 1 + e; 2) — 1 + e > w > — 1 — e; and 3) w < — 1 — e. Working in the conformal Newtonian gauge, the perturbations 
of DE can be described by 

5 = -(l + w)(0-3<f>)-3H(c 2 s -w)5 , (3) 

9 = -H(l-3w)9 —e + k 2 (^- + ^) . (4) 

1 + w 1 + w 

Neglecting the entropy perturbation, for the regions 1) and 3), the EOS does not across —1 and the perturbation 
is well defined by solving Eqs. (|3l4[) . For the case 2), the perturbation of energy density 8 and divergence of velocity, 
9, and the derivatives of 5 and 9 are finite and continuous for the realistic quintom DE models. However for the 
perturbations of the parameterizations, there is clearly a divergence. In our study for such a regime, we match the 
perturbations in region 2) to the regions 1) and 3) at the boundary and set 

5 = , 9 = 0. (5) 

In our numerical calculations we limit the range to be \Aw = e\ < 10~ 5 and find our method to be a very good 
approximation to the multi- field quintom. More detailed treatments can be found in Ref. [22| . 

The publicly available Markov Chain Monte Carlo (MCMC) package CosmoMC[58( is employed in our global fitting, 
and modifications have been made to include DE perturbations, and to suit the DE models which we study. We assume 
purely adiabatic initial conditions and work in the flat universe with VL to tai = 1- Our most general parameter space 
is: 

P = (cjb, u c , © s , r, w , w a , n s , ln(10 10 'A s )) (6) 
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Table I. Mean la constrains on the EOS of DE and f2 m . The left columns are obtained with 
"WMAP3+SDSS+2dFGRS+SN so i,;+GRBs" combinations and the right columns are correspondingly from without GRBs. 





WMAP3+LSS+SN aoW +GRBs 


WMAP3+LSS+SN 9oW 




ACDM constant w dynamical w(a) 


ACDM constant w dynamical w(a) 


Wa 


-1 -0.853^ -l.OOSt^ 
0.533+°^ 
2qn+ ' 020 285+ 021 2Q2 + ° : ° 21 

U.ZJU_ Q 020 U.Z,OO_ 2 U.z»z_ 021 


-1 -0.863ir ^ -1-OOli^ 
0.443l°;i|5 

n OCR+0-020 n noo+0.020 n 9S o+0.021 



where cut, = il^h 2 and uj c = tt c h 2 are the baryon and cold dark matter densities relative to the critical density, 
Q s is the ratio (multiplied by 100) of the sound horizon at decoupling to the angular diameter distance to the last 
scattering surface, r is the optical depth due to re-ionization, Wq and w a is the parameters of the EOS of DE, A s and 
n s characterize the power spectrum of primordial scalar perturbations. For the ACDM, wq = —1, w a = 0, 

We vary the above parameters and fit to the observational data with the MCMC method. For the pivot of the 
primordial spectrum we set A; s o = O.OSMpc" 1 . The following weak priors are taken: r < 0.8, 0.5 < n s < 1.5, 
—3 < wo < 3, and —5 < w a < 5. We impose a tophat prior on the cosmic age as 10 Gyr < to < 20 Gyr. Furthermore, 
we make use of the HST measurement of the Hubble parameter Hq = lOOh km s _1 Mpc^ 1 by multiplying the 
likelihood by a Gaussian likelihood function centered around h = 0.72 with a standard deviation a = 0.08 59]. We 
also adopt a Gaussian prior on the baryonic density VL^h 2 = 0.022 ± 0.002 (1 a) from Big Bang Nucleosynthesis 60] . 

In our calculations, we take the total likelihood to be the products of the separate likelihoods (Ci) of CMB, LSS, 
SNIa and GRB. Defining \\ $ = — 21og£i, we then have 

2 2 , 2 I 2 , 2 

Xh,total — Xl,CMB ' Xh,LSS Xh,SNIa "+" Xl,GRBs ■ (' ) 

If the likelihood function is Gaussian, x\ coincides with the usual definition of y 2 up to an additive constant corre- 
sponding to the logarithm of the normalization factor of C. For CMB, we include the three-year WMAP (WMAP3) 
data and compute the likelihood with the routine supplied by the WMAP team @. For LSS, we use the 3D power 
spectrum of galaxies from SDSS [6l[ and from 2dFGRS[62j. To minimize the nonlinear effects, we restrict ourselves 
only to the first 14 bins when using the SDSS results[63j], the range of k is 0.01578 < k/h < 0.10037. For SNe la, we 
mainly present the results with the recently released "gold" set of 182 supernovae published by Riess et al. in Ref.[3]. 
For the GRB data, we have considered the published sample by Schaefer[56j], which includes 69 GRB events. Upon 
using these distance modulus, we notice the circulation problem associated with GRBs for cosmological probing. Due 
to the lack of the local calibration, usually the correlations depend on the cosmological parameters that we attempt 
to constrain. However, there are results that show the relation does not change dramatically in a wide range of 
cosmological parameters (64j. In this paper we do not adopt the correction for the circulation problem, and we take 
the distance modulus published by Schaefer. In the calculation of the likelihood from SNe la and GRBs data, we 
marginalize over the nuisance parameter [65l [66| . 

For each regular calculation, we run 8 independent chains comprising of 150, 000—300, 000 chain elements, and spend 
thousands of CPU hours on a supercomputer. The average acceptance rate is about 40%. We test the convergence 
of the chains by Gelman and Rubin criteria[67j and find that R — 1 is on the order of 0.01, which is much more 
conservative than the recommended value R — 1 < 0.1. 



III. RESULTS AND DISCUSSIONS 

In this section, we present our global fitting results. We summarize the la constrains on the corresponding param- 
eters in Table I. We focus on the EOS of DE and £l m . In order to show explicitly the effects of GRBs, we compare the 
results between the two cases with or without GRBs. From the Table, we can see that the best fit values and the errors 
change when the GRBs data are included, especially on the parameters of the EOS of DE. For ACDM or constant 
w DE models, the best fit values and the errors change little when considering the GRBs data, however, for the 
dynamical DE models, the effects from GRBs can be obviously seen. The best fit value of (wo,w a ) is (—1.001,0.443) 
for the combined SN la + WMAP3 + SDSS + 2dFGRS. When we adopt the GRBs data, the best fit value of (w , w a ) 
changed to (—1.005,0.533), and the error bars shrink considerably, especially on w a which delineate the evolution of 
DE. One can find that the 2a error bar of w a changed from 0.443^ ( j ) '502 *° 0.553^^ which is tightened about 15%. 
This can be seen in Fig 1 graphically. 

In Figure [TJ we delineate the two dimensional posterior constraint on wq — w a . The solid lines and the dashed lines 
stand for la and 2a constraints respectively. We divide the parameter space into four regions representing different 
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dark energy models by two lines: wq = — 1 and wq +w a = — 1. For models located in the upper left region labeled as 
Quintom A, the equation of state of dark energy crosses —1 from upside down, i.e., w is greater than —1 in the past 
and becomes less than —1 at present. The evolution of w for models of Quintom B has an opposite direction. From 
the plot, it is noted that Quintom A models are mildly favored by current observational data. The AC DM model 
continues to be a consistent one at 2a level with and without GRBs included. 

The parameter w a represents the evolution of DE. It is known that the CMB data alone cannot constrain well the 
dynamics of DE. LSS data are valuable in DE studies mostly because they provide a tight limit on f2 m , which in 
turn helps to constrain the properties of DE due to the degeneracy between Vt m and the EOS of DE in cosmological 
observables. Currently the measurements of the luminosity-distance from SNe la at various redshifts give rise to the 
most direct constraints on the dynamics of DE. In a flat universe with the EOS of DE given by Eq. (1), the luminosity 
distance can be written as 



H [n m (i + z>Y + (i - n m )(i + z ')m+-o+v Ja ) exp [_ 3Wa _!_]^ ' 

From this equation, one can see that there are degeneracies between the background parameters wq, w ai Q m and Hq. 
In figure [2] we present the degeneracy by exploring the information of luminorsity distance at different redshifts in 
the wq — w a parameter plane. Different colored bands describe the parameter space of wo,w a where the variation 
of cLl is in between ±1% for z = 0.1 (black), 0.5(red),l(green), 2(blue),3(cyan),10(magenta) and llOO(yellow). One 
can find that the degeneracy between wo and w a varies with the redshift, which in turn implicit that combining 
the information of at different redshifts can help break such a degeneracy. Therefore, in order to constrain the 
cosmological parameters well, one needs distance determinations for a wide range of redshifts. For current SNe la data, 
the redshift range is limited with the highest observed redshift ~ 1.7. On the other hand, for the GRB sample used in 
our analysis, the redshift extends to as high as ~ 6.3 with 39 data points having z > 1.5, thus, the complementarity 
of GRBs to SNe la is highly expected. This plot is the ideal case for showing the degeneracy between the parameters 
wo and w a for different redshift, because we fix the other parameters except wq and w a . In fact, if we vary the other 
parameters like Q m or Hq, the changing of the degeneracy with redshift will also exist but will not be so obvious as 
figure the different colored bands will be much broader than the ones in figure [2] respectively. Figure [1] is more 
general results which is obtained from the global fitting where we have all the parameters in equation ([6]) varying. 

To demonstrate the degeneracies between different parameters, we show in Figure[3]the two-dimensional correlation 
contours of Wq, w a and Q m from our global fitting. It is seen that, with the GRBs data, the constraints on the 
parameters are tightened to a certain degree, especially for wo and w a . Our results show that high-redshift cosmological 
probes can play roles in the study of the dynamics of dark energy, and thus it is of importance to explore high-redshift 
cosmological probes. 

GRBs can bridge up the gap between the relatively nearby SN la and the much earlier CMB. Comparing with SN 
la, GRBs have their own advantages. They are very powerful and thus can be detected out to very high redshifts. 
Furthermore, the gamma-ray photons suffer from no dust extinction when they propagate to us. Therefore they 
deserve detailed studies. Although the cosmological applications of GRBs are currently limited by the relatively small 
number of available data and the quality of the data, it is considerably important to investigate their potentials and 
related problems. With the accumulation of observational data and the advances of theoretical understandings of 
GRB physics, the goal of using GRBs as cosmological candles could be better achieved. 

IV. SUMMARY 

In this paper, we have made the first global fitting with the combined GRB, WMAP3, SDSS, 2dFGRS and SN 
la data. We concentrate on dynamical DE models, and explore the complementarity of GRBs to other cosmological 
probes. DE perturbations are treated carefully. From the global fitting results given by MCMC, we find that high- 
redshift GRBs may have effects in constraining the EOS of DE, especially for the dynamical DE models. Including 
GRBs data can shrink the contours and modify the best fit values of DE parameters. The ACDM model is consistent 
with the current data within 2a confidence level, with Quintom A-like models mildly favored. 

GRBs are the most powerful astrophysical events in the universe, which hold great potential to reveal the high 
redshift universe. Even though, the GRB data are not as good as SN la currently, our results indicate the possible 
potentials of GRBs in dark energy studies. 

It is worth mentioning that our global fitting method is different from using the CMB shift parameter, the linear 
growth factor or the BAO parameter measured from the LSS. Although these parameters carry important information 
of dark energy, their specific values are usually extracted from observational data under certain conditions. Overlooking 
these conditions leads to inappropriate use of the values of these parameters, which in turn could result biased 
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FIG. 1: 2a constrains in the w and w a plane, the dotted red line is given by using WMAP3+182 "Gold" SNIa+LSS+GRBs 
with our 8-parameter parametrization discussed in the text. Solid black curves come from without GRBs. For both cases we 
considered perturbed dark energy. 




FIG. 2: The different color region simble ±1% variation around lines of constant cLl at redshift 0.1 (black), 0.5(red),l(green), 
2(blue),3(cyan),10(magenta),1100(yellow), taking A CDM model as fiducial model. This plot delineate the degeneracy between 
the parameters wo and w a at different redshift z. 




FIG. 3: 1-d posterior constraints and 2-d joint 68% and 95% confidence regions for the parameters wo, w a and Q, m obtained via 
MCMC methods. The dotted red line is given by using WMAP3+182 "Gold" SNIa+LSS+GRBs while the solid black curves 
come from without GRBs. For both cases we considered perturbed dark energy. 
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constraints on cosmological parameters. Our MCMC analysis are performed on observational data directly and 
therefore avoid such problems. We also take into account dark energy perturbations, whose effects are not included 
in the shift parameter and the BAO parameter. 

The analyse presented in this paper mainly aim at emphasizing the importance of high-redshift cosmological probes, 
which are not limited to GRBs. Our results demonstrate their contributions to dark energy studies clearly, especially 
on the dynamics of the dark energy component. 
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